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f~^ , To go beyond Gaussian approximation to the Hohenberg-Kohn free energy playing the key role 

f^ ' in the density functional theory (DFT), the density functional integral representation would be 

C~^ \ relevant, because field theoretical approach to perturbative calculations becomes available. Then 

CN| . the present letter first derives the associated Hamiltonian of density functional, explicitly including 

5_^ ' logarithmic entropy term, from the grand partition function expressed by configurational integrals. 

CL(, Moreover, two things are done so that the efficiency of the obtained form may be revealed: to 

•^r ' demonstrate that this representation facilitates the field theoretical treatment of the perturbative 
calculation, and further to compare our perturbative formulation with that of the DFT. 

^^ ! PACS numbers; 64.10. +h, 83.70.Hq, 05.40. +j, 05.70.Ce 

The Hohenberg-Kohn (HK) free energy — in terms of the density functional theory (DFT) — is the natural exten- 
^ . sion of the Helniholtz one. At the beginning, we see what is meant by this via the primary definition of the HK free 
c/3 ; energy [||. 

Let us start with a grand canonical system which has a volume of V and is surrounded by a reservoir of a chemical 
potential fi in the unit of the thermal energy fc^T ||l| . We consider one-component classical particles, to keep notations 
as possible as simple, and define the grand potential U, in the form: 
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with r being the position vectors of particles, J the external potential, and U{ri, rj) the two-body interaction potential 

^^ for particles i and j. In the thermodynamic limit, the Helmholtz free energy F of a canonical system with N particles 

JZ^ [ can be obtained from 51 by a Legendre transform that F = il + /i7V, where N is the averaged total number given by the 

—^ ■ relation N — —dfl/d^ in the grand canonical formalism. Similarly, the HK free energy Fhk is defined by replacing 

(^ . N in this mapping with the averaged density field ip{r) given in terms of the shifted external field J{r) = J{r) — /i as 

o 

a^ ■ <P = (p) - -~, (2) 

a^■ ^■^ 



where p = X]i=i '^('" ^ ^0 i^ the density operator and (• • •) represents the ensemble average with the weight in Eq. 
(nl) : the Legendre transform of fi with use of ip and J yields 



"O ■ Fhk{<p)^W)-V-J. (3) 
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where ip ■ J = drip J. Since the relation (ph is reciprocal, Fhk exactly satisfies the equation that SFhk / Sip = —J. 

In spite of the above generality and tractability of the HK free energy Fhk, the introduction of Fhk has been 

mostly in the context of the DFT. Indeed the DFT is one of the most powerful tools for investigating not only spatially 

inhomogeneous states for simple liquids but also a variety of more complex fluids (e.g. liquid crystals and polymers) 

j^ [01 . However, it is also to be noted that there are cases where the largely fluctuating systems such as fluids near 

critical points are beyond its scope. This is seen from the following Ramakrishnan-Youssouff form [p|: 



AFi^= fdr^^{r)\n^-Udrdr'C(^\r-r';pl)Aip^{r)A^^{r'), 
J Pm ^J 



(4) 
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where A F^j^ is the excess free energy around an arbitrary uniform density p in an A-domain (e.g., Hquid region in 
hquid-vapor coexisting state), C''^'>{r — r';p ) is the second-order direct correlation fmiction, and Aip^ — ^p ~'p is 
the density difference between p and Lp obtained from the self-consistent equation: 

<^(r)=exp[C(i)(r;^)- J(r)]. (5) 

The merit of this form (0) is that short scale correlation — particularly crucial for fluids — may be taken into elaborate 
account via the input of the direct correlation function; this insert is much of benefit due to the extensive study on 
the Orstein-Zcrnike integral equation B. For all that, the limitation is also to be realized that the direct correlation 
functions as input effectively consider only quadratic fluctuations in case these are the solutions of the integral equation 
within the mean spherical approximation including the equivalence, i.e., Percus-Yevick one to hard-core potentials [pi. 
A systematic way of going beyond the Gaussian approximaion to the HK free energy is to start with the density 
functional integral representation so that field theoretical approach to the perturbative calculations may be available. 
Then we purpose, first of all, to transform straightforwardly the primary definition of Fhk expressed by configurational 
integrals to the functional integral representation: (i) to derive from Eq. (|^) with the relation (^ the following 
expression that 

e- ^"" ^ I Dp exp [-{ H,ad{p) - V ■ J}] (6) 

Hsadif)^ drdr' -p{r)U{r,r')p{r')+p{r)J + p{r)lnp{r)-p{r), (7) 

where p denotes an instantaneous quantity (not operator) of density. 

The justification of the Hamiltonian (|^) is the main result of this letter; one has had few grounds so far for the form 
([7|), except for the primitive discussion that division of ideal gas systems into cells produces the entropy term, i.e., 
the last two terms on the right hand side of Eq. (m) MM , though the above expression M) comprising the familiar 
free energy functional is trivial intuitively and hence the corresponding form for the Helmholtz free energy has been 
often used a priori ^^ ■ 

(i) 

As usual, let us insert into the position vectors' form (IT]) of the grand partition function, the identity with use of 
the auxiliary field ilj{r): J dpdil) exp [i'0(p ~ p) ] = 1- Then Eq. (||) reads 

e"^^ = I Dp LIV exp [ ~H (p, ?/') ] (8) 

H{p,^) — drdr'-pUp' + pJ — ip'ip~exp{—i'iJj + p), (9) 

with setting p{r) — p and p{r') = p'. Previous treatments of Eqs. (||, |^) have conventionally proceeded to perform 
Gaussian integration over p for fixed ip — the Hubbard-Stratonovich transformation, what is called |l^,|l^. The 
reduction is exact indeed, and hence there appears to be no choice but to do so. As a matter of fact, however, 
Gaussian approximation — this also considering quadratic contribution — would be applicable to iJj for given p: we have 
the alternatives in mapping Eqs. (PP) to more tractable forms. 
Then, taking the latter approach, Eqs. (|8|j|) are reduced to 

e-'' = J Dpexp[~H,ad] (10) 

and Hsad{p) given by Eq. (|^); Hsad is equal to the Hamiltonian (Q) along the saddle point path for tp, i.e., H{p,ipsad) 
with ijjsad satisfying the saddle point equation that 5 H{p,^j))/54'\'4)=Tp,^a — 0- To be noted in Eq. (pX|), the excess 
grand potential Af7 = — In [J Z?^ ^-pi^i') /2] arising from quadratic fluctuations oi 6ip — ip — tpsad is absent. This is 
due to the following evaluation: Gaussian integration over ^ — carried out by discretized flelds, pi and ipi — yields the 
apparently nontrivial term of the Lee- Yang type [|2| such that 

An= lim ^ fdrUn{pia^), (11) 



in the continuum limit (or the vanishing hmit of the lattice constant a defined as a^ — (po)~^ = '^/N), where the other 
trivial term has been formally absorbed into the integral measure Dp following the standard procedure |12| . The excess 
potential Afl given by Eq. ([ll|), however, converges to 0.5(A'' — N), the half of the difference between the actual total 
number N and the most probable value N, as found from the expansion that \n{pi a^) = {pi ~pa) a^ + 0[{{pi —pa) a'^}^]- 
Thus AD, is negligible in the thermodynamic limit. 

In the next step, we equate the functional derivative 6il/6J using the representation ([1C|) with the averaged density 



ip obtained from both the relation (0) and the position vectors' expression (m) of fJ: we put ip = (p) — {p)c with (• ■■)c 
denoting the average under the weight in Eq. (|lO|). Then the Legendre transform (^) leads to the expressions (Jg, 0), 
in question, of the HK free energy Fhk', the principal purpose of this letter has been accomplished. 

What we have to do, in addition, would be to show virtues of the now justified representations (|g, [7|). Then two 
things are done in the remainder: (ii) to demonstrate that this form facilitates the field theoretical treatment of the 
perturbative calculation, and (iii) to compare our perturbative formulation with that of the DFT. 

(ii) 

Let us return to the grand potential il given by Eq. (UG). The saddle point equation that 5Hsad{p)/ ^ p\p=p . j=o = 
in the absence of external potential, J{r) = 0, produces the mean-field density: 

p^ (r) = exp I - y dr' U{r,r')p^ {r') + p^. (12) 

Expanding around p the logarithmic term in the Hamiltonian difference AH = Hsad{p) — Hgadifj^,)-, we obtain 

^-ci^^-h^Mpm) f Dp exp{~{AH + p-J)} (13) 



AH{p)^ fdrdr'lp{U+^-^^^^)p' + Eip}, 
J 2 p^ 



(14) 



where p = p ~ p , the fluctuation of the total number is neglected (i.e., / drp « 0), and -E(p) denotes the terms 
higher than quadratic ones due to the logarithmic expansion. 

The Hamiltonian AH given by Eq. (|l^ has some characteristics other than standard field theoretical formulations 
[ p2[ : One is that the free part of AH includes the positionally dependent coefRcient 1/p (r), indicating that our 
formalism is applicable to investigating density fluctuations in structured fluids where the mean-field density itself is 
not uniform but spatially oscillates. Next, AH reveals that higher terms E(j)) in one-component systems arise from 
entropy allowing thermally activated hopping processes and not from interactions pq ]. Finally it is noted that J is 
reduced to J. This is why the external potential J is not included in determining the mean-field density from the 
saddle point equation; if we take as a reference density the mean-field one in the presence of J, AH has no external 
source to create a generating functional. 

Of particular interest is the first prop erty; this, for example, makes it possible to generallize Debye-Hiickel equation 
as Fisher et al. propose recently ||l^,|lj]. The present letter, however, restricts itself to the simplified systems which 
have only the short-ranged potentials and consist of macroscopic subsystems (or domains) inside which the mean- 
field density p (r) take constant values except the boundaries. A typical situation is the coexisting phase of two 
subsystems (such as liquid and gas). 

Consequently, neglect of the interfacial energy reduce AH to the sum of the contributions from the domains. With 
the volume and the mean field density in an s-domain (s = A, B, • • ■) denoted by Vg and p" , respectively (where the 
overbar is for emphasizing the constancy) , the functional integral in Eq. (|l3) reads in terms of the Fourier transformed 
density, p'^{k) = (l/Vg) J dr exp{ik ■ r)p''(r): 
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(15) 



where further shifted density is abbreviated to p* , C is a constant independent of J, and the propagator A|^ is given 
in the usual form (A|^)^^ = wq + (1/p^ ) + W2fc^ due to the conventional expansion of the short-ranged potential, 

U = uq + U2k^ , in the fe-space. 



The expression (|_5|) implies that Feynman graphs are now available. Thus the HK free energy Fhk gets into 
spotlight, because Fhk given by Eqs. (||J3) includes the generating functional of the ra-point vertex functional F^"^ 
consisting only of one-particle irreducible diagrams as well as the standard field theory |1^: putting that Aip'^ = ^^Ptf 
in an s-domain {s ^ A, B , ■ ■ •) as in Eq. (Q) and ignoring the number fluctuation (i.e., / dr Aip'^ « 0) as before, Fhk 



Fhk^ E Pmf + ^F^k (16) 

fdrdr'-p' Up''+-p' Inp" -p' (17) 

J dr, ■ --dr,, r(")(ri, • ■ • ,r„; r^) A^^(ri), • ■ • , A^^(r„). (18) 
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It is to be noted in the above representations that Fhk correctly includes the mean-field Helmholtz free energy Fj^ip 
in the absence of external potential J, in contrast to previous density functional integral formulations ||ll[ which take 
as a reference the free energy for the density smeared over the entire system. As a result, the present formalism may 
focus on appropriate density deviation from not a smeared value but a mean-field one p* , and therefore makes the 
perturbative approximation more precise than previous approaches. 

(iii) 

Even at the starting point, there exists difference between the DFT and our formalism: the former separates the 
entropy term, J dr iplncp — ip, from Fhk beforehand, whereas the latter does the mean-field free energy Hsadip^r) 
(see Eq. (p^)). Such distinct types of theories, though, have correspondence in some cases as will be seen below. 

Most successful is the absence of interactions between particles, i.e., C/ = 0, where the field theoretical representation 
is identified with the DFT formalism. To see this, we consider the grand potential fl. In the DFT, fi is reduced to 
n = J dr (fiax^iKfiex + 'fiexJ with tlic avcragcd density (pex = exp(— J) obtained from the relation (|5|) when C^^' — 0. 
While the saddle point equation, SHsad/Sp = 0, explicitly including the external potential J also produces the density 
equal to i^cx, and therefore the deviation of Q in the functional integral formalism from that for the DFT is given in 
the first evaluation as — In [/ D Sp exp ( —(5pY /2lPcx) ], being the same kind as the excess grand potential AfJ. Then 
repeating the similar discussion to that after Eq. (|ll|), the difference may be ignored and thus the equivalence in the 
case of f/ = is assured. 

For L/ ^ 0, on the other hand, let us compare both representations of the excess HK free energy AF^^ in an 
A-domain up to quadratic terms for Atp^. For the DFT formulation (13), the expansion of the logarithmic term yields 

AF^^=y'drdr'l('-C(2) + fc^'| A^^{r)A^^{r'). (19) 

In applying the mean spherical approximation (MSA) to the calculation of the direct correlation function for hard- 
core(- Yukawa) or square well fluids |j], we have only to put — C^^-* = [/ in Eq. ([l9|), because the other condition that 
the pair distribution function is to be set to zero inside hard spheres with the diameter d is formally satisfied due to 
the potential of C/ = cxd for r — r' < d. On the other hand, the functional integral representation ( p^ ) under the tree 
approximation reduces to AF^j^ — AH^Aip"^) |lj], and therefore takes the same form up to quadratic terms as the 



DFT expression (19) with use of the MSA for the above mentioned fluids (see also Eq. (|l4|)). 

The conformity for [/ 7^ contrarily results in highlighting some merits of the field theoretical forms (|l^) to (|lq), 
which we describe in conclusion. One virtue other than the DFT using the MSA is that the density functional integral 
representation may still take fluctuations into more elaborate consideration systematically, by including loop graphs 
in Eq. (Iq). Moreover, we would like to stress that a reference density — this being merely an arbitrary value in the 
DFT (see the statement just after Eq. (Q)) — is obtained from the relation ( p^ ) self-consistently in our formalism. 
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